<N 



> 
>< 



Non-linear noise excitation 
and intermittency under high disorder^ 

Davar Khoshnevisan Kunwoo Kim 



Abstract 



o 

(N 

University of Utah University of Utah 

c3 

^ Last update: March 3, 2013 

in 

4^ Consider the semilinear heat equation dtu = d^u + Xa{u)^ on the 

^ interval [0 , 1] with Dirichlet zero boundary condition and a nice non- 

^ random initial function, where the forcing ^ is space-time white noise 

'— ' and A > denotes the level of the noise. We show that, when the 

solution is intermittent [that is, when inf^ |cr(2:)/z| > 0], the expected 
^ L^-energy of the solution grows at least as cxp{cA^} and at most as 

exp{cA^} as A — > oo. In the case that the Dirichlet boundary condition 
is replaced by a Neumann boundary condition, we prove that the i^- 
energy of the solution is in fact of sharp exponential order exp{cA^}. 
We show also that, for a large family of one-dimensional randomly- 
forced wave equations, the energy of the solution grows as exp{cA} as 
A — )• 00. Thus, we observe the surprising result that the stochastic 
wave equation is. quite typically, significantly less noise-excitable than 
its parabolic counterparts. 
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1 Introduction 



The principal aim of this paper is to study the non-hnear effect of noise in 
stochastic partial differential equations whose solutions are "intermittent." 
Our analysis is motivated in part by a series of computer simulations that 
we would like to present first. 

Consider the following stochastic heat equation on the interval [0 , 1] with 
homogeneous Dirichlet boundary conditions: 



lUt{x) = -^:^ut{x) + \ut{x)^ for < X < 1 and t > 0, 



d , , 1 

dt ^ ' 2dx^ 
ut{0) = ut{l) =0 for t > 0, 

uo{x) = sin(7rx) for < x < 1. 



(1.1) 



Here, ^ denotes space-time white noise on (0,oo) x [0,1], and A > is an 
arbitrary parameter that is known as the level of the noise. 

The stochastic partial differential equation (1.1) and its variations are 
sometimes referred to as parabolic Anderson models. Those are a family of 
noise-perturbed partial differential equations that arise in a diverse number 
of scientific disciplines. For a representative sample see Balazs et al [1], 
Bertini and Cancrini [2] , Carmona and Molchanov [8] , Corwin [9] , Cranston 
and Molchanov [10], Cranston, Mountford, and Shiga [11,12], Gartner and 
Konig [17], den Hollander [15], Kardar [20], Kardar et al [21], Majda [23], 
Molchanov [24], and Zeldovich et al [28], together with their substantial 
bibliographies. 

Since the initial function is the principle eigenfunction of the Dirichlet 
Laplacian on [0,1], it is easy to see that the solution is exactly ut{x) = 
sin(7r,T) cxp{— 7r^t/2} when A = 0; this is the noise-free case, and a numerical 
solution is shown is in Figure 1. 

Figures 2, 3, and 4 contain simulations of the solution to (1.1) for respec- 
tive noise levels A = 0.1, A = 2, and A = 5. 

One might notice that, as A increases, the simulated solution rapidly 
develops tall peaks that are distributed over relatively-small "islands." This 
general phenomenon is called intermittency; see Gibbons and Titi [18] for 
a modern general discussion of intermittency and its ramifications, as well 
as connections to other important topics [in particular, to turbulence]. Our 
discussion is motivated in part by the discussion in Bliimich [3] on the quite- 
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1 ^. 




Figure 1: No noise (A = 0.0, maximum height = 1.00) 



1.4 ^. 
1.2 ^. 
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Figure 2: Small noise (A = 0.1, maximum height w 1.05) 
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Figure 3: Modest noise (A = 2.0, maximum height ~ 30.53) 
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different topic of NMR spectroscopy. 

Presently, we study non-linear noise-excitation phenomena for stochastic 
partial differential equations that include the parabolic Anderson model (1.1) 
as a principle example. Namely, let us choose and fix a length scale L > 0, 
and then consider the stochastic partial differential equation, 



d 



dt 



ut{x) 



- ^ 

ut{L) = 



ut{x) + Xa{ut{x))£, for < X < L and t > 0, 

for t > 0, 



(1.2) 



where a : R — t- R is a Lipschitz-continuous function, ^ denotes space-time 
white noise as before, and the initial function uq : [0 ,L] R+ is a non- 
random bounded continuous function that is non-negative everywhere on 
[0 , L] and strictly positive on a set of positive Lebesgue measure in (0 , L). 
For the sake of simplicity, we assume further that a{0) = 0, though some of 
our work remains valid when cr(0) 7^ as well. 

It is well known that the stochastic heat equation (1.2) has an a.s.-unique 
continuous solution that has the property that 

sup sup E (\utix)\''^ < 00 for all T > and k e [2 , 00). (1.3) 

xelO,L] t6[0,T] ^ ^ 

We will be interested in the effect of the level A of the noise on the [expected] 
energy S't{\) of the solution at time t; the latter quantity is defined as 



#,(A):= Je(k||2,j^^^,) (t>0). 



Throughout, we use the following notation: 



inf 
2eR\{0} 



Lo- := sup 
zeR\{o} 



cr z 



(1.4) 



(1.5) 



Clearly, ^ £a ^ Lo-- Moreover, L^- < 00 because a is Lipschitz continuous. 

The following two theorems contain quantitative descriptions of the non- 
linear noise excitability of (1.2). These are the main findings of this paper. 
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Theorem 1.1. For all t> 0, 



^ ^ liminf log4(A), limsup ^ log4(A) ^ 8L^t. (1.6) 

In a companion paper [22] wc show that cxpIcA"*} is a typical lower 
bound for the energy of a large number of intermittent complex systems. 
In particular, the energy of the solution for the stochastic heat equation on 
[0 , L] with a periodic boundary condition is shown to be of sharp exponential 
order exp{cA^}. The following theorem says that we can get the same kind 
of result when wc replace a Dirichlet with a Neumann boundary condition, 
provided additionally that the initial profile remains bounded uniformly away 
from zero. 

Theorem 1.2. Suppose that we replace the Dirichlet boundary condition in 
(1.2) by a Neumann boundary condition; that is, we suppose that u solves 
the following stochastic heat equation: 

d d"^ 

— -nt(x) = Tz-^ut(x) + Xa(ut(x))^ for < x < L and t > 0, 
ot ax^ (17) 

If, in addition, we assume that iTiixG[o,L] uo{x) > 0, then for every t > 0, 

f^t 1 1 Ql/f 

^liminf -log 4(A) ^limsup -log 4(A) ^ -f^. (1.8) 
o7re AtcxD A^ _\-|-Qo A^ 16 

Theorems 1.1 and 1.2 together show that under fairly natural regularity 

conditions [that include £a- > 0], the energy behaves roughly as cxp{cA^}, 
which is a fastly-growing function of the level A of the noise. In Section 4 
we document the somewhat surprising fact that, by contrast, the stochastic 
wave equation has typically an energy that grows merely as exp{cA}. In 
other words, the stochastic wave equation is typically substantially less noise 
excitable than the stochastic heat equation. 
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2 Proof of Theorem 1.1 



As is customary, we begin by writing the solution to the stochastic heat 
equation (1.2) in integral form [also known as the mild form], 

ut{x) = {PtUQ){x) + X I pt-s{x,y)a{us{y))^{dsdy), (2.1) 

J{0,t)x{0,L) 

where {Pt}t^o denotes the semigroup of the Dirichlet Laplacian on [0,L] 
and {pt}t>o denotes the corresponding heat kernel. That is, in particular, 
Poh = h for every h e L°°[0 ,L], and 

{Pth){x) := [ pt{x , y)h{y) dy for alH > and x£[0,L]. (2.2) 
Jo 

One can expand the heat kernel pt{x,y) — ^the fundamental solution to the 
Dirichlet Laplacian on [0 , L] — in terms of the eigenfunctions of the Dirichlet 
Laplacian as follows: 

oo 

Pt{x ,y):=Yl (kn{x)4>n{y)e~^-'\ (2.3) 

n=l 

where 

1 /2 

/Xn := j , Mx) ■■= [^ij sm (— j , (2.4) 

for all n ^ 1 and ^ x ^ L. According to the maximum principle: (i) 
Pt{x , y) > for alH > and x,y E (0 , L); and (ii) pt{x ,y)dx < 1 for all 
y G [0 , L] and t > 0. 

Next we recall, briefly, how to one establishes the existence of a mild 
solution to (1.2). 

Let uf'\x) := uq{x) for all x G [0,L] and t ^ 0, and then define itera- 
tively, for all A; ^ 0, 

ui''+'\x) = {PtUo){x) + X [ pt-s{x,y)a{ui''\y))^{dsdy). (2.5) 

J{0,t)x{0,L) 

Then the method of Walsh [27, Chapter 3] — see Dalang [13] — shows that 
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{u^"^}t^o is locally uniformly Cauchy in i^^(P) in the sense that for all T > 0, 



V sup sup Je(\4''^'^\x)-uI''\x)\A <oo. (2.6) 
f.^ixe[o,L]te[o,T] V ^ ^ 

It follows fairly easily from this that u := liink-^^ u[''\x) exists in I/^(P) 
and solves (2.5). Moreover, one can deduce that u is unique among all mild 
solutions that satisfy (1.3) — see Dalang [13] for the details. 

Proof of the first bound in Theorem 1.1. We follow an idea that is classical 
in the context of PDEs (see Kaplan [19]); in the context of SPDEs, a non- 
trivial adaptation of this idea was used in Bonder and Groisman [4] in order 
to compare the solution of an SPDE to the solution to a one-dimensional 
diffusion in order to describe a certain blow-up phenomenon. In a somewhat 
loose sense, we follow the same general outline, but need to make a number 
of modifications along the way. 

Recall that denote the eigenfunctions of the Dirichlet Laplacian 

on [0,L] and {^iJ,n}'^=i are the corresponding eigenvalues. Since pt{x ,y) = 
Pt{y ,x), an appeal to a stochastic Fubini theorem [27, Theorem 2.6, p. 296] 
implies the following for all n ^ 1 and t > 0: With probability one, 

{ut , (f>n) 

= {uo , PtM + A / (Pt_s0„)(y)(T(«s(y)) ^(ds dy) 

J(0,t)x(0,L) l^-'J 

= e-'^"*(uo , 0„) + A / e-'^«(*-^)0„(y)(7(u,(y)) ^(ds dy). 

J(0,t)x(0,L) 

Consequently, the Walsh isometry [27, Theorem 2.5, p. 295] shows us that 

E[(ut,</.„)2] =e-2/^"*(no,</.n)' 

+ A2 fds f dye-2'^"(^-^)|</>„(y)|2E(|a(u,(2/))p) 
Jo Jo V y 

+ \Hl j\-^^^-^'-'Us dy \cl>r.{y)\^E{\us{yf) • 
We may apply the Cauchy-Schwarz inequality, together with an appeal to 
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the Fubini theorem, in order to see that 



^ |0„(t/)|2E(|u,(y)|2)dy5J ^E[K,<^„)2]. (2.9) 

Thus we see that, for every fixed n ^ 1, the function 

F(t):=e2^"*E[(^xt,0„)2] (t > 0) (2.10) 
satisfies the recursion 

F{t) ^ {uo , (l)nf + X^il [ F{s) ds for all t > 0. (2.11) 

Jo 

Thus the Gronwall's inequality [or just recursion, in this case] shows that 
F{t) ^ {uo , (pnf exp {X^ilt) (t > 0). (2.12) 

Equivalently, 

E [{Ut , (Pnf] > {uo , <Pn? exp ( [X^l - 2fXn] t) {t>0). (2.13) 

Since X]^i('Ut , 0n)^ ^ ll^*lli2[o l]' thanks to Bessel's inequality, it remains 
to prove that the following is strictly positive: 



:= ^{uo , (pnfe-^^^^ for all t ^ 0. (2.14) 



n=l 



But Q > for all t ^ simply because we can write ct = \\PtUo\\'j^2^Q j^y and 
this expression is strictly positive, thanks to the maximum principle and the 
fact that «o > on a set of positive measure. This concludes the proof of 
the lower bound on the energy in Theorem 1.1. □ 

Our derivation of the upper bound of Theorem 1.1 requires a few ele- 
mentary preparatory lemmas. 

Lemma 2.1. For all /3,t> 0: 

sup /* e-^(*-^) ds dy \pt-s{x , y)]^ ^ (2.15) 
xe[o,L]Jo Jo ^vP 
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and 

rt 



^ e-^(*-^)ds^ dxj^ dy \pt-s{x,y)? ^ (2.16) 
Proof. Because \(f)n{x)\ ^ (2/L)V2 for all x G [0 , L] and n ^ 1, 

ft f-L ft t» 

/ e-'^(*-)d. / dy\pt-six,y)f= J2i^n{x)fe-'^f'^^-^^'-'^ ds 
Jo Jo Jo 

1 °° 1 

n=l \ / / 

The first bound follows because 

^ 1 f"" da; _ L 

for all /3 > 0. In order to obtain the second bound, we integrate (2.17) [dx], 
using the fact that ||0n||L2[o,L] = 1- D 

We are now ready to establish the energy upper bound in Theorem 1.1. 

Proof of the second hound in Theorem 1.1. We appeal to the nonlinear renewal- 
theoretic technique of Foondun and Khoshnevisan [16] in order to establish 
the following a priori estimate: For every 5 G (0 , 1), 

E(|u,(x)r) ^ r™/2||«o||roc[o,Liexp (2tm' (r=^)') ' ^^"^^^ 

valid uniformly for all m G [2,oo), x G [0,-L], and t ^ 0. Once this is 
established, then we set m = 2 in (2.19), and then integrate [dx], in order to 
deduce the upper bound of the theorem. 

Since L^- < oo and a{0) = 0, it follows that 

\aiz)f ^ hl\z\^ for all z G [0 , L]. (2.20) 

Recall that u"^^ denotes the fcth-stage Picard iteration approximation to u. 
According to the Burkholder-Davis-Gundy inequality (see Burkholder [5], 
Burkholder et al [6], and Burkholder and Gundy [7]; see also Foondun and 
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Khoshnevisan [16] for an explanation of these particular constants), 



(0,t)x(0,L) 



^ 4m 



Jo Jo 



dy \\Zs 



(2.21) 



for all predictable random fields Z, all t > 0, and m e [2 , oo), where ||Q||m := 
{E(|Qp)}V"^ for all random variables Q. Therefore, we can deduce from 
(2.5) that 




^ II^oIIl-io.l] + Ay4mL^ / ds j dy [ptsix ,y)]'^ uy{y) 
Define 



iVfe(;3):=sup sup (e"'^* u'j.''\x) 

t^O x£[0,L] 



1/2 



The preceding and Lemma 2.1 together show that 



ds I dy[pt-s{x,y)f uf\y) ^[N^m 



l2 6' 



fit 



Therefore, we combine our efforts in order to deduce the bound 

uf^^\x) ^ \\uo\\L^yoM + V2mAWe^*/2^-V4iVfe(^). 



(2.22) 



(A; ^ 0, /3 > 0). (2.23) 



(2.24) 



(2.25) 



Since the right-hand side is independent of x G [0 , L] and depends on t 
only through the term exp(/3t/2), we can divide both sides by exp{pt/2) and 
optimize over t and x in order to deduce the recursive inequality. 



(2.26) 
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So far, /3 was an auxiliary positive parameter. We now fix it at the value 



:= (2.27) 



where J G (0 , 1) is arbitrary but fixed. For this particular choice, 

Nk+iiP,) ^ ||no||Lco[o,L] + (1 - 5)NM {k ^ 0). (2.28) 
Because Nq(P^) = ||'Uo||l°°[o,l]! preceding iteration yields 

oo 

supiVjk(^*) ^ ||«o||l'x>[o,l] y](l - Sy = 5"^||iio||L°°[o,L]- (2-29) 

This and Fatou's lemma together imply that 

sup sup (e~'^'^\\utix)\\lA S'^\\uo\\Looio n, (2.30) 
t^o xe[o,L] ^ ^ 

which, in turn, yields (2.19). □ 

3 Proof of Theorem 1.2 

Let us now consider the following stochastic heat equation on the interval 
[0 , L], subject to a Neumann boundary condition: 

d 92 

—ut(x) = —-^uAx) + Xaiutix))^, for < X < L and t > 0, 

Ot OX'^ (n 1) 

where ^ denotes space-time white noise as before, and the assumptions on a 
and uq{x) are also the same as for the Dirichlet problem (1.2). For the sake 
of simplicity, we assume further that cr(0) = 0, though some of our work 
remains valid when (7(0) 7^ as well. We may assume that 

4 > 0. (3.2) 
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Otherwise, the more interesting part of Theorem 1.2, that is the lower bound 
on the energy, has no content. [The proof of the upper bound will not require 
(3.2).] Also, let us define 

£ := inf Unix). (3.3) 
xe[o,L] 

According to the hypotheses of Theorem 1.2, e > 0. 

As wc did for the Dirichlet problem, we begin by writing the solution to 
the stochastic heat equation (3.1) in integral form. Namely, we write 

ut{x) = {PtUo){x) + X Pt-six , y)a{us{y)) ^(ds dy), (3.4) 

J{0,t)x{0,L) 

where {Pt}t^o denotes the semigroup of the Neumann Laplacian on [0 , L] 
and {pt}t>o denotes the corresponding heat kernel. That is, Poh = h for 
every h E L°°[0 ,L], and 

{Pth){x) := [ pt{x , y)h{y) dy for alH > and xe[0,L]. (3.5) 
Jo 

It is well known that we can write the heat kernel for {Pt}t^o in the following 
form: 

oo 

Pt{x,y):= [rt{x-y-2nL)+rt{x + y-2nL)], (3.6) 

n=— oo 

where F denotes the fundamental solution to the heat equation in R; that 
is, 

1 / z^' 



^'^^^''^7^t'''''^[~Tt' (*>0'^eR). (3.7) 



[This requires an application of the method of images.] 

The existence and uniqueness of a solution of (3.1) is shown in Walsh 
[27, Chapter 3]. The following estabhshes the lower bound of Theorem 1.2. 

Proposition 3.1 (An energy inequality). For all t > 0, there exists a finite 
and positive constant K — independent of X — such that 

<ft(A)^lexpf%^V (3.8) 



K " V 87re 

simultaneously for every A ^ 1. 
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Proof of Proposition 3.1. Owing to (3.6), the Neumann heat kernel is conser- 
vative; that is, Jq Pt{x ,y)dy = 1. Therefore, we apply the Walsh isometry 
to find that for all fixed e,t > and x G [0 , L], 

^(\ut{x)\^)^e^ + \^ j\s j\y \pt-s{x,y)f¥.{\a{us{y))\^) 
^ + {Ma? j\s j\y \pt-s{x,y)fY.{\us{yf) . 

The preceding is a recursive, nonlinear renewal-type, inequality for the 
function (t , x) i— >■ E(|'Ut(x)p). Even though this renewal inequation is the 
lower bound's counterpart to (4.11), we cannot solve this inequation using 
the renewal-theoretic ideas of Foondun and Khoshnevisan [16]. The reason 
for this is that the method of [16] works for large values of the time variable 
t\ whereas here we have t fixed, but A large. Instead, we appeal to old ideas 
of localization of multiple integrals [i.e., classical large deviations] that are 
probably due to P.-S. Laplace, though we could not find explicit references. 
Namely, we expand the self-referential inequality (4.17) and observe that 
most of the contribution to the resulting multiple integrals occur in a small 
portion of the space of integration. 

Now we carry out our program by first expanding, to one term, our 
renewal inequality as follows: 

e(|«,(x)|2) 

^ £2 + e''{\iaf f ds [ dy \pt-s{x , y)? (3.9) 
Jo Jo 

+ (A4)^ fds f dy \pt-s{x,y)? f dr [ dz \ps-r{y , (Wr{^)f) 
Jo Jo Jo Jo ^ ^ 

and proceed. In order to see better what is happening, let us introduce a 
family of linear operators V as follows: For all t > 0, a; G [0 , L], and all 
Borel- measurable functions ijj : {0 , oo) x [0 , L] — >■ R+, 

, x) := /* ds r dy \pt-s{x , y)Ms , y). (3.10) 
Jo Jo 

Also define l(i , x) := 1 for alH > and x G [0 , L]. Then our recursion can 
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be written compactly as follows: 

oo 

E(|«t(x)|^) ^£2^(A4)''=(P'=l)(t,x); (3.11) 



k=0 



where P° := 1 and := V'^V for all A; ^ 0. We integrate (3.11) [dx] in 

order to obtain the following key lower bound: 

oo „i 

l^tWf^e^yZiXiaf' iV''l){t,x)dx. (3.12) 

We examine the preceding sum, term by term, using induction. 
The first term in the sum is trivial; viz., 

/ iV°l){t,x)dx = L. (3.13) 
Jo 

Therefore, let us consider, as a first step, the second term in the sum, which 
is straight forward: 

/ iPH)it,x)dx= [ ds [ dx 1 dy\pt-s{x,y)f -.= ^{1). (3.14) 
Jo Jo Jo Jo 

Before we estimate the third term of the infinite sum in (3.12), let us observe 
that 



L 

2 



{V^l){t,x)dx (3.15) 

/^t /"L /'L f'S f'L 

= ds dx dy \pt-s{x,y)f dr dz \ps-r{y,z)f. 
Jo Jo Jo Jo Jo 
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Therefore, we can bound that third term, from below, as follows: 



/ {V'^l){t,x)dx (3.16) 
Jo 

ft f'L f'L f'S f'L 

^ ds dx dy \pt-six,y)f dr dz \ps-riy,z)f 

Jt/2 Jo Jo Jo Jo 

^ ds dx dy \pt-s{x ,y)]'^ / dr dz \ps-r{y , z)]"^ 

Jt/2 Jo Jo Js-t/2 Jo 

pt pL pL r^/*^ 
^ ds dx dy [pt-six,y)f / dr dz \priy,z)f 

Jt/2 Jo Jo Jo Jo 

t-t/2 i-t/2 i-L i-L pL 

= ds dr dx dy \ps{x ,y)f / dz \priy,z)f. 
Jo Jo Jo Jo Jo 



It is not important that we have a t/2 and in the bounds of our integrals on 
the second line; t/e would have worked equally well for us. The key point 
is that most of the contribution to our multiple integral comes about where 
s ~ t and r ~ s in the second line. 

Now, let us recall that pr{y , z) can be thought of as the transition func- 
tion for Brownian motion in (0 , L) , reflected upon reaching the boundary 
{0 , L} [26] . The Markov property of reflected Brownian motion implies 
the semigroup property of {Pt}t^o and in particular that fQ[pr{y , z)]"^ dz = 
P2riy ,y)- And by symmetry, jQ\ps{x,y)]'^dx = P2s{y ,y), as well. Conse- 
quently, we obtain 



/ iVH){t,x)dx^ [ 
Jo Jo 



t/2 



\p2siy,y)fds 



dy. 



(3.17) 



Once we understand this, we can see how to bound the remaining terms in 
the infinite sum in (3.12) as well. In fact, induction shows that 



fl. pt/{k+l) 

/ iP''+'l){t,x)dx^ / / \p2siy,y)?ds 
Jo Jo Jo 



k+1 



dy, (3.18) 



for all A; ^ 1. Once again we stress that our multiple integrals are bounded 
from below by a very small portion of the region of integration. 
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Finally, we apply Jensen's inequality in order to see that 



Jo Jo Jo 



k+1 



(3.19) 



for all A; ^ 1. Recall (3.14) in order to see that 



Jo 



^ {t/k) 
L 



(3.20) 



for all A; ^ 1 and t > 0. This and (3.12) together yield 



k=i 



L 



(3.21) 



Finally, we estimate the function $ from below. 
For all T > 0, 



f'T f'L f'L 

$(t) = / ds / dx- / dy \ps{x ,y)f 
Jo Jo JO 

= ds dx p2s{x,x). 
Jo Jo 



(3.22) 



Since P2s{x , x) ^ r2s(0) = (Stts) for all s > — see (3.6) — the preceding 
yields the pointwise bound 



(3.23) 



Thus, we obtain the following: 

oo 



k=l 

oo 



oo 

j = l L 




(3.24) 



47re 
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thanks to the elementary bound, {j/ey ^ j\, vaHd for all integers j ^ 1. 
The proposition follows from the preceding and the Taylor series expansion 
of the exponential function. □ 

Before we prove the upper bound in Theorem 1.2, let us state a simple 
result regarding the Neumann heat kernel. 

Lemma 3.2. For every e > there exists a positive and finite constant 
K := K^^L such that 

rt 3 + e 

sup sup / e~^^p2s{x ,x)ds^ for all (3^ K. (3.25) 

Proof. In accord with (3.6), for every s > and x G [0 , L], 

oo 

P2s{x,x)= [T2s{'^nL) + V2s{2x-2nL)] 



^3r2,(0) + 2£]r2,(2nL)+ r2,(2|nL - a;|) 



n=l neZ: 

\nL—x\^l 



(3.26) 



^3r2.(0) + C(L) 
= ^ + C{L), 

uniformly for all s > and x £ [0 ,L], where C (L) is a positive and finite 
constant that depends only on L. The numerical bound of 3 [in front of 
r2s(0)] accounts for the fact that, depending on the value of x, there can 
be at most two choices of n G Z such that \nL — x\ < 1. We integrate the 
preceding in order to find that 

e-^'p2s{x ,x)ds^-^ + ^, (3.27) 



which has the desired effect. □ 

Wc need only to prove the upper bound in Theorem 1.2; the correspond- 
ing lower bound has already been established. 

Proof of the second bound in Theorem 1.2. According to the Walsh isome- 



18 



try, for alH > 

E{\ut{x)\') (3.28) 

= \{PtUo){x)f + X'' dsj^ dy\pt-six,y)]^E(^\a{uM)f) 

^ lho|lioc[o,L] + (AW)2^ ds dy \pt-s{x,y)]^E(^\usiy)\^^ . 

We solve this inequality by applying the method of Foondun and Khosh- 
nevisan [16]; namely, let us define, for all /? > 0, 



^(^):=supsup e-f^^E {\ut{x)f) 



(3.29) 



The preceding and Lemma 3.2 together imply that for all £ > there exists 
K := K^^L such that for all P ^ K, 

=^(-5) ^ ll^olliooro n +^(/3)(AW)2sup sup / e-^'\p2s{x , x)f ds 

(3 + e)(AW)2 



We let < (5 < 1 and choose /3 as 



(3.30) 



(3 + £)^(AW)^ 
8(1 - (5)2 ' 



in order to deduce the inequality 



(3.31) 



(3.32) 



valid as long as A is large enough to ensure that /3* ^ K. In this way we 
find that 



E{\ut{x)f) ^ ;f||^^o||ioc[o,L]exp 



i3 + efi\K)H 
8(1-5)2 



(3.33) 



We integrate both sides of (3.33) [dx] in order to obtain the upper bound, 
after some relabeling. □ 
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4 A stochastic wave equation 



In this section we consider the non-hnear stochastic wave equation 

-^wt{x) = -^wt{x) + Xa{wt{x))^ for X G R and i > 0, (4.1) 

subject to non-random initial function wo{x) = and non-random non- 
negative initial velocity vq G L^{^) H -^^(R) such that ||fo||L2(R) > 0. It is 
well known (see Dalang [13] for the general theory, as well as Dalang and 
Mueller [14], for the existence of an L^-valued solution) that there exists a 
unique continuous solution w to (4.1) that satisfies the moment conditions, 

sup supE (\wt{x)\'') + sup S't{X)<oo, (4.2) 

for every G [2 , oo) and T > 0. 

The main result of this section arc the following bounds on the energy 
^^(A) := {E(||ut||2,^j^p}i/2 of the solution to (4.1). 

Theorem 4.1. For every t > 0, 

^ liminf ilog4(A) ^ limsup^log^t(A) ^ (4.3) 



4v8 -^^oo A A^oo A V8 

Define 

2 

air) ^= ^ /\ J > 0)- (4-4) 

The preceding theorem requires the following elementary real- variable fact. 

Lemma 4.2. Suppose g is a non-negative element o/L-'^(R)nL2(R). Then, 
there exist positive and finite constants Ai and A2 — depending only on ||5'||li(r) 
and llffllia^R) — such that 

A^H{t) ^ f dz f dy {g * g) {y - z) ^ A2H{t), (4.5) 

J-t J-t 

for all t > 0, where g{x) := g{—x) for all x G R. 

Proof. Define h := g * g for simplicity, and note that: (i) h ^ 0; (ii) 
h G L^(R)nL2(R) with ||/i||li(r) = ll5'llii(R) and h{0) = ||5'|||2(r)- Further- 
more, since h(x) is maximized at a; = 0, thanks to well-known facts about 
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continuous, positive-definite functions. Now we put these facts together in 
order to see that 

f dy [ dzh{y-z)^4h{0)t^A2t\\h\\Li(^) 
J-t J-t (4.6) 

^ 16 A H{t). 

This proves the upper bound [with an exphcit ^2]- 

On the other hand, J^^dy J^^dz h{y - z) = {4 + o{l))t^h{0) as t 4, 
and dy /^^ dz h{y — z) = {2 + o(l))i||/i||^i(R) as i t 00. The lower bound 
follows from these observations, since h is non-negative. □ 

Proof of Theorem 4-1- The solution to the stochastic wave equation (4.1) 
can be written in mild form as follows: 

wt{x) = \Wt{x) + \X f l[o,t-s\{\x-y\)<y{ws{y))adsdy), (4.7) 

J(0,t)x'R. 

where ^ 

Wt{x) := j VQ{x-y)dy. (4.8) 

Therefore, the Walsh isometry for stochastic integrals assures us that 
E(K(x)|2) 

nt poo (4 9) 

= l\Wtix)f + lX'' ds dyl[o,t-s]{\x-y\)E{\a{ws{y))f), 
Jo J-00 

whence by Fubini's theorem, 

ft f'OO 

Wt{X)f = l\\Wt\\hiK) + h>^' / (i-s)dW dyE{\a{ws{y))f) 

Jo J-00 

f't poo ^ ' 

^ill^tlli2(R) + 5^'L' j/^-'^"^' J_^dyE{\ws{y)\'). 

Since ||W''t|||2(R) = jltdy j!i^dz {vo*vo){y — z), Lemma 4.2 ensures that 
the squared energy |<ft(A)|^ of the solution to the stochastic wave equation 
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satisfies the renewal inequality, 

|^t(A)|' ^ ^2i'||^^o|li2(R) + IX'^Ll j\t - s) {Ss{X)f ds, 



for all t > 0. 
Since 



the preceding implies that 

4^2 



sup 



(e/3)2' 



(4.11) 



(4.12) 



^ (^ll^o|li2[o,i] + iA2L2^(/3)^*(i-s)e-^(*-^)d. 



|2 



(4.13) 



2/32 



Let us choose and fix an arbitrary S G (0, 1). We can apply an a priori 
method in order to show that < 00, where 



/3* := 



(4.14) 



See the derivation of (2.19); we omit the details. In this way, we can solve 
(4.13) in order to deduce the following: 



8^2||t'o|li2[o^^] 

^(eAW)2 



(4.15) 



Equivalently, 



^ ^(eAL.)2 ^^^iTiiraj 



(4.16) 



This readily yields the upper bound of the theorem. 

We prove the corresponding lower bound by observing the following coun- 
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terpart of (4.11), which holds for the same reasons as (4.11) does: 



\^tiX)\'^A,H{t)\\vo\\h[o,L] + '2>^'^' 



\t-s)\<^sWf ds. (4.17) 







Now we proceed as we did in the proof of Proposition 3.1; we expand our 
renewal inequation as an inequality in terms of an infinite sum of multiple 
integrals of increasingly-high powers. And then show that the multiple in- 
tegrals are large, when A S> 1, mainly because of the contribution of the 
integrand in a small region of integration. In this way the remainder of the 
proof is exactly the same as the derivation of Proposition 3.1. However, as 
it turns out, one has to be quite careful in order to guess the correct region 
of integration, as it will be very significantly larger than the one in the proof 
Proposition 3.1. 

With the preceding in mind, we begin by writing 



where y := (t — si) x (si — S2) x • • • x — Sn) X H{t), over the range of 

the integral in (4.19). Let us emphasize that the nth term involves an n-fold 
integral that is integrate on a large part of the original region of integration; 
this is in sharp contrast with the proof of Proposition 3.1, where the effective 
region of integration was extremely small for the nth term when n S> 1. Once 
the correct region of integration is identified, we can continue the argument 
in the proof of Proposition 3.1. Namely, we obtain the following bounds. 




etc. In this way, we obtain the following generous lower bound. 



l<^*(A)P 



(4.19) 
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after we appeal to time reversal: 

oo 

|4(A)r^Ai||t;o||i.[o,L]^^W-E(§^'^') 

/.si/2 |.s„-i/2 n {4:.20) 

/ dsi / ds2 • • • / ds„ TT 

JO JO Jo 



X / dsi / dS2 • • • / ds„ H Sj. 



A change of variables and induction together yield 

H/2 i-si/2 /-s„-i/2 n 



/ dsi / ds2--- dsn T\sj 

Jo Jo Jo 

ft rsi fSn-l " 

= 4-" / dsi / ds2 • • • / ds„ n s,- (4.21) 
Jo Jo Jo 

~ 4" • (2n)!!' 

where (2n)!! := (2n) x (2n — 2) x • • • x 2 denotes the usual double factorial 
of 2n. We require only the simple bound (2ra)!! ^ (2n)! in order to deduce 
that 

|4(A)r^^i||.o|li.[o,Li^(*)-E(2^J (2rz)!- ^'-''^ 

Since 

for A ^ 2\/8/(tO. It follows that 

|4(A)|^ ^ i^i||^o|li.[o,^]i?W • E (^J (4-24) 
whenever A ^ 2y/8/{ticr)- Because 

E(^)'^^=0(A) as a too, (4.25) 
the lower bound follows. □ 
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